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Several origins have been proposed for the Equatorial Highlands on Venus, includ- 
ing (i) spreading centers, analogous to mid-ocean ridges on Earth (Head and Crumpler, 
1987), and (ii) plume-related uplift, analogous to hotspots on Earth (Phillips and Ma- 
lin, 1984; Stofan and Saunders, 1990). Recently, the spreading center hypothesis has 
been shown to be incompatible with the measured geoid and topography variations 
over the Highlands (Kiefer, 1990). It is also difficult to reconcile the range of geoid 
anomalies over the Highlands with a steady-state plume model. There is a large varia- 
tion in admittance values (geoid/ topography ratios) among Highland regions (see table 
1). This variation suggests different uplifted regions represent distinct stages in a 
time-dependent process. Herrick and Phillips (1990) have proposed that the Highland 
regions listed in table 1 are formed by large mantle diapirs. According to this model, 
topography and geoid height decrease with increasing age of the highland, as the diapir 
spreads out beneath the lithosphere. A similar dynamical process has been proposed for 
continental flood basalt and ocean plateau volcanism on Earth (Richards et al, 1989). 
On Earth, rapid uplift and massive volcanism are hypothesized to accompany forma- 
tion of a new hotspot, when a large mantle plume diapir rises toward the lithosphere 
and undergoes decompression partial melting. 

In order to determine if the diapir model is compatible with the sequence of tectonic 
and volcanic events recorded in the surface geology of the Highlands, a series of finite 
difference calculations have been made of the ascent and partial melting of a spherical 
thermal diapir in am incompressible, infinite Prandtl number, isoviscous fluid. The 
upper surface of the fluid is a rigid lid, representing the crust. The initial conditions 
include a conductive thermal boundary layer at the upper surface, representing the 
thermal lithosphere, and a buoyant spherical diapir with initial temperature Tp, located 
near the base of the fluid. The diapir is released at the beginning of the calculation, 
and allowed to rise. The diapir shown in Figure 1 began with an excess temperature 
AT = Td — T oo =400K. Appropriate thermodynamic and transport parameters were 
chosen for the mantle of Venus. Partial melting is calculated using a simplified pressure- 
dependent melting law. The dynamical behavior of the melt phase is represented as 
a bimodal process. Where the porosity <f> is below a critical value, <f c , percolation 
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Table 1: Equatorial Highlands 


Dome 

Diameter (km) 

Elevation (km) 

Geoid (m) 

Beta Regio 

2800 

4.5 

90 

Thetis Regio 

3000 

4 

70 

Ovda Regio 

3400 

4.8 

35 

Artemis Plateau 

2500 

1.5 

<20 


Table 2: Diapir tectonics 


Stage 

Time (Myr) 

Tectonic Style 

Volcanic Style 

1 

0-5 

Rapid Uplift + Extension 

Rift Volcanism 

2 

5-20 

Rapid Subsidence 

Flood Volcanism 

3 

>20 

Slow Subsidence + Marginal Compression 

none 


Figure 1: A numerical calculation of the evolution of a rising thermal diapir. Upper 
panels: Profiles of geoid height (solid) topography (dotted) and dynamical to- 
pography (dashed). Lower panels: Contours of the temperature field. Shaded 
regions show partial melt concentration. 
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